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Phoretic self-propulsion: a mesoscopic description of reaction dynam- 
ics that powers motion 1 
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The fabrication of synthetic self-propelled particles and the experimental investigations of their dynamics have stimulated interest 
in self-generated phoretic effects that propel nano- and micron- scale objects. Theoretical modeling of these phenomena is often 
based on a continuum description of the solvent for different phoretic propulsion mechanisms, including, self-electrophoresis, 
self-diffusiophoresis and self-thermophoresis. The work in this paper considers various types of catalytic chemical reaction at the 
motor surface and in the bulk fluid that come into play in mesoscopic descriptions of the dynamics. The formulation is illustrated 
by developing the mesoscopic reaction dynamics for exothermic and dissociation reactions that are used to power motor motion. 
The results of simulations of the self-propelled dynamics of composite Janus particles by these mechanisms are presented. 



1 Introduction 

Nano- and micron- scale motors that are able to propel them- 
selves through solution by utilizing the energy derived from 
chemical reactions constitute an interesting class of objects 
with numerous potential applications. Such small self- 
propelled objects have been made in the laboratory and their 
properties have been investigated. In particular, a good deal 
of research has been devoted to some of the earliest exam- 
ples of motors of this type, bimetallic and striped metallic 
rod motors 1-6 . Experimental studies of Janus colloidal mo- 
tors, where one face of a polystyrene or silica sphere is coated 
with platinum 7 ' 8 , and sphere dimer motors, where two linked 
spheres, one silica and the other platinum 9 , have also been 
carried out. Interest in such motors stems from their potential 
applications 4 ' 10-12 . For example, striped metallic rod motors 
that include Ni can be controlled by external magnetic fields 
and experiments have demonstrated that targeted cargo deliv- 
ery can be achieved 13-15 , as has transport by Pt-Au nanomo- 
tors using photochemical stimuli 16 . Self-propelled Janus par- 
ticles have recently been investigated experimentally for cargo 
transport applications 17 ' 18 . 

Such motors are placed in a solution containing fuel of some 
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type, most often hydrogen peroxide in the cases mentioned 
above, and chemical reactions on portions of the motor give 
rise to directed motion. In these examples the motors move 
by self-phoresis, where the gradient of some field across the 
motor, which is generated by asymmetrical chemical reactiv- 
ity, induces fluid flow in the surrounding medium resulting 
in propulsion. The metallic rod motors exploit an electro- 
chemical mechanism involving oxidation and reduction reac- 
tions at the rod components, leading to self-electrophoresis 5 , 
while the Janus and sphere dimer motors move by self- 
diffusiophoresis. Motors that are propelled by the generation 
of bubbles have also been constructed and studied 19-21 . While 
many aspects of propulsion by self-phoresis are generic, the 
motor motion depends on the nature of the chemical reactions 
taking place on the motor, as well as those in the environment 
in which the motor moves. 

A characteristic feature of such self-propelled motion is that 
the entire system is force-free. This has implications for mo- 
tor propulsion and for the nature of the flows that develop in 
the system. The treatments of self-propulsion based on the 
motion of colloidal particles in gradients rely on various con- 
tinuum approximations 22-26 . While these approximations are 
certainly valid for large particles, and even for quite small par- 
ticles, as the nano and molecular levels are approached such 
approximations should be subject to scrutiny. There have 
been numerous studies of propulsion by self-phoretic mech- 
anisms 27-32 . In common with our earlier studies of sphere 
dimer and polymer motors where a simple reaction A — >> B oc- 
curs on the catalytic sphere 33-38 , we consider a particle-based 
mesoscopic description of the dynamics of self-propelled par- 
ticles. The dynamics is such that the mass, momentum and en- 
ergy conservation laws are satisfied and hydrodynamic flows 
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are correctly described. 

In this article we show how the mesoscopic dynamics can 
be extended to treat a general class of chemical reactions on 
the catalytic sites of the motor. The reactive dynamics is con- 
structed to satisfy the conservation laws so the basic prop- 
erties of the dynamical laws are preserved. As examples of 
these more general reactions, we present the results of simu- 
lations of self-propulsion of a composite Janus particle aris- 
ing from a dissociation reaction A — »• 2B and an exothermic 
reaction A — > B + Au at motor catalytic sites. These exam- 
ples illustrate some interesting effects; in particular, in the lat- 
ter case of an exothermic reaction propulsion occurs by self- 
thermophoresis. The general strategy described in this paper 
should provide a framework that is useful in explorations of 
the effects of complex chemical kinetics on self-propelled mo- 
tion at a mesoscopic level. 

The outline of the paper is as follows: The system we 
consider comprises a motor and the environment in which it 
moves. In Sec. 2 we describe the motor itself, the (possibly 
reactive) environment in which it moves and how it interacts 
with the chemical species in the environment. The detailed de- 
scription of the reactive processes at the motor catalytic sites is 
given in Sec. 3 and the results of simulations of the motor dy- 
namics are presented in Sec. 4. The conclusions of the paper 
are given in Sec. 5. 

2 Mesoscopic dynamical system 

The chemically-powered motor we study is a composite Janus 
particle comprising a rigid assembly of catalytic and noncat- 
alytic beads (see left panel of Fig. 1). One half of the Janus 
particle is made of catalytic (C) beads while the other is half 
is made from non-active (N) beads. The solvent in which the 
Janus particle moves contains, in general, reactive and nonre- 
active chemical species. The beads in the Janus particle in- 
teract with solvent particles through specific intermolecular 
forces. There are no explicit solvent- solvent intermolecular 
forces. These are treated by multiparticle collisions 39 ' 40 as 
discussed below and described in reviews 41 ' 42 . In addition, the 
catalytic beads are able to catalyze various chemical reactions 
that are responsible for the propulsion of the motor. It is the 
nature of these chemical reactions that is the principal focus of 
this investigation and the description of how general reactions 
are implemented on the catalytic sites will be the topic of the 
next section. Here we specify the nature of the mesoscopic 
dynamics that governs the evolution of the system. 

The evolution of the entire system, the Janus particle and 
the surrounding solvent is carried out by combining molec- 
ular dynamics (MD) with multiparticle collision dynamics 
(MPCD) 39 ' 40 . In this hybrid evolution scheme, Newton's 
equations of motion are integrated for all interacting particles 
for time intervals T. This includes the Janus particle and all 
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Fig. 1 Left panel: A representation of the bead model for the Janus 
particle. The top part is the catalytic side C of the Janus particle and 
the bottom part is the noncatalytic side N. Right panel: The 
coordinates system defined with the Janus particle's center of mass 
as the center. The gray disk represents the Janus particle. In 
subsequent figures, that region contains empty data. 

solvent particles that interact with it. The remainder of the 
noninteracting solvent particles free stream during these time 
intervals. At the time intervals % (possibly) reactive multipar- 
ticle collisions take place as described in the next subsection. 

2.1 Dynamical model for the environment in which the 
motor moves 

The motor environment is modeled by reactive multiparticle 
collision dynamics (RMPCD) 43 , since we wish to allow for 
the possibility that chemical reactions may occur in the bulk 
of the solution 36 in addition to the catalytic reactions that 
occur on the motor itself. RMPCD combines multiparticle 
velocity-changing collisions among solvent species with local 
probabilistic rules for changing the identities of those solvent 
species that participate in reactive events in the bulk of the 
solution. To carry out the dynamics, at the time intervals 1 
the system is partitioned into cells labeled by £, . We suppose 
that various reactions, specified by the index j, occur in the 
solution: 

Lv£x a ^£v£x« (7 = l,...,r). (1) 

The stoichiometric coefficients for reaction j are and 
while kj and k-j are the rate constants for the forward and 

reverse reactions. Each reaction, £ a V^X a -4 L a v^X a , is 

taken to occur with probability pj (N^ ) = dj(l — e~ a o T ^j ja^ , 

where = (n\ ,Af| , • • • , ) is the set of numbers of the dif- 
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ferent species in cell and aj = kj(V c )Y[ a N aU ( N ot ~ v a)- 
accounts for the number of different ways the reaction can oc- 
cur in the cell, with = L/tfj?. The rate constants kj(V c ) 
are scaled to account for the cell volume V c . Reaction rules 
may be constructed so that mass, momentum and energy are 
conserved in the reactive events 43 . 

Once the local reactive events have taken place, all particles 
in each cell £ undergo multiparticle collisions where the post- 
collision velocity of particle i in cell v-, is given by 39 

▼5 = V 5 + fi> 5 -(v/-V 5 ) (2) 

where is the center of mass velocity of all particles in cell 
£ and co^ is a rotation operator for the cell £ . This algorithm 
preserves mass, linear momentum and energy within each cell. 
The combined effect of multiparticle collisions and reactive 
events is to mimic the effects of many real reactive and non- 
reactive collisions that take place in the time interval T. Since 
the full hybrid MD-RMPCD preserves the basic conservation 
laws, hydrodynamic interactions and hydrodynamic flows in- 
duced by the self-propulsion are properly described so that 
this mesoscopic simulation scheme can capture all essential 
features of the Janus particle and solvent flow motions. 

3 Chemical reactions at motor catalytic sites 

In this section we show how different types of reaction can 
be implemented at a mesoscopic level at the catalytic face of 
the Janus particle. As noted above, in order to respect the ba- 
sic conservation laws and insure that fluid flows are properly 
described, the reactions must be constructed to satisfy mass, 
momentum and energy conservation laws. In earlier studies 
on sphere dimer motors we restricted our considerations to the 
simple thermoneutral catalytic reaction, A + C ^ B + C. This 
reaction amounts to a "coloring" process that changes A to B 
on encounters with the catalytic C sphere. Mass, momentum 
and energy are simply conserved in this process. The propul- 
sion arose from the fact that the A and B species interacted 
with different intermolecular potentials at the noncatalytic N 
sphere. 

The actual catalytic reactions that take place on motor ac- 
tive regions may be quite complex, involving adsorption onto 
the surface with subsequent formation of products and their 
release, or binding to active sites in enzymes followed by con- 
formational or other changes and product release. All of these 
events could involve dissipation of energy into local solvent 
degrees of freedom, either directly or subsequent to various 
internal energy transfer processes within the motor. Our meso- 
scopic treatment of such reactive events coarse grains over 
most of these molecular-level details. In particular, when the 
reagents lie within the interfacial zone around the catalytic 



portion of the Janus particle where intermolecular forces be- 
tween the reagents and Janus particle beads are non-zero, there 
is a possibility that a reaction will take place. We subsume all 
of the reaction details into effective reactive events where re- 
active species conversion and possible energy distribution in 
solvent degrees of freedom take place just as the species exit 
the reactive interfacial zone. The local solvent molecules that 
participate in any energy exchanges with reactive species are 
taken to lie in the MPC collision cell that contains the reactive 
species. Using this strategy it is possible to construct a general 
class of catalytic reactive events on the active face of the Janus 
particle 44 . 

In this coarse-grain description of real reactions catalyzed 
on the Janus motor, the conservation of momentum and energy 
in a general reactive event can be written as 

m JVj + £ mvi = mjv'j + £ mv'i, 0) 

for the momentum, where £ denotes the cell in which the re- 
acting particles are found and the prime again denotes post- 
reaction quantities. Summing over all particles in the cell £ 
allows one to utilize the velocities of solvent particles, other 
than the reacting ones, to satisfy conservation rules. Energy 
conservation reads 

+ \ L m i y2 i + L u i = \ m J yf J + \ L m ^'i + L u 'f 

(4) 

As described below, in order to model endothermic or exother- 
mic reactions we may associate an internal degree of freedom 
with the chemical species and in this equation u a is the in- 
ternal energy of species a. For thermoneutral reactions this 
internal state label is not needed. 

Below we shall make use of the following notation: m a is 
the mass of species a; M a p £ is the sum of the masses of the 
different species, for instance, M a p = m a + mp and M aa = 
2m a \ r a is the position of the particle a; R a p is the center of 
mass position of particles a and j3 ; r a p = (r a — rp )/ |r a — | 
is the unit vector along the a-f5 direction; \ a is the velocity 
of the particle a; \ a p = (m a \ a + mp\p)/M a p is the center 
of mass velocity of a and j3 ; w aj g = \ a — is the relative 
velocity of a and j3; ji a p = m a mp/M a p is a reduced mass; 
and u a is the internal energy associated to a particle a. For 
some of the quantities, for instance the sum of the masses or 
the center of mass velocity, arbitrary arrangements of indices 
can be made. For quantities where the order matters, such as 
w aj s or jA a p, grouping can be made with a parentheses; e.g., 
Va(py) ™ a mp r /M a py. 

The precise manner in which reaction rules are constructed 
depends on the chemical mechanism one wishes to simu- 
late. We illustrate the procedure by considering two differ- 
ent reaction schemes: an exothermic unimolecular reaction 
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A + J ^ B + J + Au, where Au is the energy release on re- 
action, and a dissociation reaction A + / — >> / + 25, where / 
represents the Janus particle. 

3.1 Exothermic unimolecular reaction 

The exothermic unimolecular reaction involves the reactant 
and product species, A and B, respectively, interacting with the 
Janus particle /. The reactive A and B species, while treated 
as structureless objects, nevertheless carry internal energy la- 
bels ua and ub, respectively. For example, this is a simpli- 
fied description of a catalyzed isomerization reaction where 
reactants and products differ by a free energy difference of 
Au = ub — ua- Mass balance requires that the masses of A and 
B be equal: ma =ntB. Recall that the Janus particle is a rigid 
object so that no Janus internal degrees of freedom need to 
be taken into account. (This condition can be relaxed when 
flexible self-propelled particles are considered.) 

There is considerable freedom in the way one chooses to 
implement a reactive scheme that accounts for the conserva- 
tion conditions in Eqs. (3) and (4). For example, for A + C — >> 
B + C + Au, we may account for the energy change Au on re- 
action by changing the post-collision relative kinetic energy 
w ksi °f me P r °duct B molecule and a solvent S\ molecule, 
and not the Janus particle. In this case the momentum conser- 
vation equations take the form, 



mj\j 
m A \A^rn Sl \s l 



mj\j 



(5) 



The second of these equations can be satisfied if the post- 
collision velocities are given by 



VA + L S 

m B 

vs x L s. 

m Sl 



(6) 



where s remains to be determined by energy conservation. 

Supposing ua>ub for now, we find that the norm of the rel- 
ative post-collision velocity W BJ depends on the pre-collision 
relative velocity waj and also on the difference of internal en- 
ergies. As v/ remains unmodified by the reaction, the depo- 
sition of the excess energy Au is made in the relative velocity 
W BS . Equation (4) becomes, in the center-of-mass velocity 
frame of B-S\, 

^bs^Is^^bs^bs, = \MBs l VBS l + \vBS{n 1 BS l +ua-u b . 

(7) 

Since the center-of-mass velocity \bs x is not altered during 
the event, Eq. (7) becomes after some algebra 



V BS l 



^ua — ub 



(8) 



From Eq. 6 we have W BSl : = ytas\ + s - Using s = |s| and 



- we find 

s 



S = 



-IWASJCOSG,- 




cos0 r ) +2 (9) 



where B r is the angle between s and v?bs x • We decide to always 
pick the + solution. The angle of s may be chosen at random. 

As a consequence of the local increase in kinetic energy an 
inhomogeneous temperature gradient will be generated in the 
vicinity of the Janus particle. If the interaction energies of the 
A and B species with the Janus particle are same, this model 
exhibits the interesting feature that propulsion will be driven 
by thermophoresis alone. 

3.2 Dissociation reaction 

Here we suppose that a dissociation reaction, A — >> B\ +B2, 
is catalyzed by interaction with the rigid Janus particle. The 
symbols B\ and B2 label the two B species molecules that re- 
sult from the dissociation of A. Mass conservation requires 
that the mass of a B molecules is half that of A: tub = rriA/2. 
Momentum conservation reads, 



mj\j + m A \ A = mj\j + m B (V Bl + \' Bl ) . 



(10) 



The energy conservation equation, Eq. (4), can be written in 
the form, 

^MA/W A/ = 2M0Bi52)/ W (£i£ 2 )-/ + 2^ BlBl ^ B ^ Bl ' ^ 

We may choose to set ^ Bi b 2 )j ~ ^' ^ e '' we ma ^ use a ^ °^ 
the relative kinetic energy of the system to change the relative 
energy of #i#2- The direction of w^ i52 may be chosen at 
random. The norm of w^ i52 is then given explicitly by 



l w flifl 2 l 



JlAJ_ 



(12) 



Other variants of these collision rules may be constructed. 

4 Simulation of Janus particle motion 

In this section we present the results of simulations of the 
self-propelled motion of the composite Janus particle when 
it catalyzes the exothermic A — >> B unimolecular reaction and 
A — >> 2B dissociation reaction. The simulations were carried 
out in a cubic box with periodic boundary conditions. This 
simulation volume contained N = Na + Nb molecules. The 
A and B particles interact with all beads in the Janus particle 
through repulsive Lennard- Jones potentials, 



V(r)=4e((-Y 



■r), (13) 
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where r c = 2 1 / 6 <7 and is the Heaviside function. Note that 
we choose the interactions of both species with both catalytic 
and noncatalytic beads to be identical; consequently the parti- 
cles are identical apart from species labels and no propulsion 
of the Janus particle can arise from the simple A — >> B ther- 
moneutral reaction considered in previous studies of sphere 
dimer motors. As described earlier, interactions between the 
A and B particles in the environment are accounted for by mul- 
tip article collisions instead of direct intermolecular forces. In 
the MD portions of the evolution, Newton's equations of mo- 
tion are integrated using the velocity Verlet algorithm with a 
time step At. The Shake and Rattle algorithms are used to 
maintain the rigid bond constraints in the composite Janus par- 
ticle. In the MPC portions of the evolution, grid shifting is 
applied at each MPCD step to insure Galilean invariance 45 . 

All quantities reported below are expressed in simulation 

units: mass ma, length a, energy £ and time y } ^f- = *sim- I n 
these units the system parameters were: simulation box size, 
48 x 48 x 48; MPC time 1 = 0.5; average number of parti- 
cles per MPC cell, 10; reduced temperature, ksT = 1/3; and 
At = 0.01 . The choice of the time step depends on the criterion 
At <C tu where tu is the natural time for the Lennard- Jones 

parameters, tu = \J The Janus particle is made from 36 
individual spheres, each with Lennard- Jones radius a = 3 and 
mass 50, so that mj — 1800. All realizations of the dynamics 
started with 100 MPCD steps with an Andersen thermostat 46 
at chosen temperature ksT = 1/3 in systems with all A parti- 
cles. We note that for these parameters the thermal velocity 

for the Janus particle is v t h = ~ 0.014. 

In order to sustain self-propelled motion in the steady state 
reactants must be supplied to the system and products re- 
moved. In the past this has been done by introducing fluxes 
of these species at the boundaries or far from the motor, or by 
bulk far-from-equilibrium reactions which are themselves as- 
sumed to be sustained by fluxes. (In simple descriptions these 
species feeds can be accounted for in pool chemical species 
that may be incorporated into effective rate constants.) In this 
study we maintain a supply of fresh reagents by taking the re- 
verse of the reaction on the motor to occur in the bulk of the 
fluid, but with different reaction rates. The reverse reaction 
for the exothermic reaction is B — >• A — Au while for the dis- 
sociation reaction it is 2B — ^ A. The reaction rates of both of 
these reactions were taken to be k r = 0.001. These reactions 
are implemented using the RMPCD algorithm 43 . Momentum 
and energy are obeyed at the cell level as in Eqs. (3) and (4) 
but without the Janus particle. Energy addition in a cell is 
obtained by a velocity scaling in the center-of-mass velocity 
frame. 



4.1 Exothermic A — » B reaction 

As noted earlier, since the A and B particles are identical apart 
from species labels, the only possible source of directed mo- 
tion (assuming instabilities do not occur) is the temperature 
gradient across the Janus particle generated by the exothermic 
reaction. In simulations starting from all A molecules in the 
simulation box, after a transient of 3000 time units, the com- 
bined effect of reactions at the Janus particle and in the bulk 
leads to a steady average value of Na and Afe. For the system 
in this steady state regime, the number concentrations ca and 
cb are displayed in Fig. 2 and show a strong excess of B at the 
top (catalytic side) of the particle. The polar axis of the Janus 
particle is used to define a coordinate system that moves along 
the particle. The definitions of the radius r and angle 6 are 
shown in Fig. 1 (right). The temperature field T(r,6) around 
the compound is shown in Fig. 3. The temperature is com- 
puted cell- wise at regular intervals during the simulations and 
binned at the center of mass of the corresponding cell. One 
clearly observes a difference between the C and N sides of the 
Janus particle, a sign that a thermal gradient is present in the 
stationary regime of the system. 




Fig. 2 The concentrations field for A and B molecules, ca and cb for 
Au = 0.5. 

The velocity of the self-propelled Janus particle can be de- 
termined from the time and ensemble average of V z (t), the in- 
stantaneous velocity of the center-of-mass velocity projected 
on the instantaneous unit vector from the center of mass TV 
hemisphere to the center of mass of the C hemisphere, 

(V z )(t) = {\£dt'r CN (t')-Vj(t')), (14) 

where the angular brackets denote an average over realizations 
of the evolution. When the time argument is omitted, (V z ) = 
(V z )(t max ) where t max is the duration of the simulation. 

Figure 4 displays (V z )(t). After a short transient, the 
running average stabilizes around an approximately constant 
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Fig. 3 The temperature field around the Janus particle for Au = 0.5. 
The black region is related to the very low solvent concentrations 
c a (r, 6) close to the Janus particle that leads to very poor statistics. 

nonzero value, indicating the existence of directed motion. 
The value of the velocity is positive, consistent with the results 
obtained by Yang and Ripoll 47 where the thermal gradient is 
generated externally on the fluid around one sphere of a sphere 
dimer. The directed velocity V z is not negligible compared to 
the thermal velocity v t h and the directed component of the mo- 
tion plays an important role. Thus, we have demonstrated that 
the chemically induced temperature gradient by the chemical 
reaction leads to self-thermophoretic propulsion of the Janus 
motor. 
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Fig. 4 The running average of the directed velocity for Au = 0.5. 
The result is obtained from an average of 16 independent 
realizations of the dynamics. 

The behavior for different values of Au is summarized in 
Table 1 where (V z ) and its standard deviation, which charac- 



terizes the coherence of V z between different runs of the same 
physical conditions, are given. As expected, increasing the 
energy Au associated to the reaction increases the propulsion 
velocity. 





Aw = 0.1 


Au 0.2 


Au = 0.3 


Au = 0.5 


(V z ) x 10 3 


1.5 


2.8 


3.8 


5.7 


a (v z ) x 1q3 


0.8 


0.9 


1.2 


1.4 



Table 1 The average directed velocity and the associated standard 
deviation for the exothermic reaction for several values of the 
exothermicity. Results are averages over 16 independent 
realizations. 



4.2 Dissociation A — » 2B reaction 

A similar set of simulations was carried out for the disso- 
ciation reaction. Figure 5 displays the concentration fields 
c a (r, 6) for the A and B species, showing the production of B 
on the top catalytic side of the Janus particle. As each A pro- 
duces 2B, the maximum value of cb is much higher than the 
average density of A particles in the system, about 13.5 par- 
ticles per MPC cell. The B particles diffuse from the C side 
of the Janus particle to the bulk of the fluid where the reverse 
reaction takes place. 




Fig. 5 The concentration fields, c«(r, 0), for A and B molecules for 
the dissociation reaction.. 

Figure 6 shows (V z )(t) [from Eq. (14)]. The gradient of ca 
and cb results in directed motion with the C hemisphere of 
the Janus particle in front. The value of V z , however, is much 
lower than the thermal velocity v t h of the Janus particles and 
the propulsive character can only be observed through aver- 
aging. Thus, even though the A and B particles interact with 
beads of the Janus particle with the same intermolecular po- 
tentials, the excess of B particles produced in the reaction pro- 
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Fig. 6 The running average of the directed velocity for the 
dissociation reaction. 

vides a mechanism for self-diffusiophoresis and propulsion is 
observed. 

5 Conclusions 

The directed motions of small nano- and micron- scale mo- 
tors that operate by chemically-fueled phoretic mechanisms 
depend on the nature of the catalytic chemical reactions that 
occur on the motor. These reactions generate the gradients of 
fields that are responsible for the propulsion of the motor. The 
work described in this investigation demonstrated how a vari- 
ety of chemical mechanisms can be incorporated into a meso- 
scopic description of the dynamics of the system. One of the 
important aspects of propulsion by phoretic mechanisms is the 
correct description of the fluid flow fields that are induced by 
the motor motion. In order to properly describe the hydrody- 
namics of the motor environment, the mesoscopic dynamics 
must preserve the basic conservation laws of mass, momen- 
tum and energy conservation. The hybrid MD-MPCD dynam- 
ics employed in this study satisfies these criteria. We have 
shown how different catalytic reactions at the motor surface 
can be constructed that maintain these conservation laws. 

The two examples we have chosen to illustrate the method 
present interesting features. For both the exothermic uni- 
molecular and dissociation reactions we have chosen to con- 
sider the situation where the A and B reactive species interact 
with the same intermolecular forces with the Janus particle. 
This is an special but interesting case that highlights the role 
that the reaction mechanism plays in the motor propulsion. 
In the case of the exothermic reaction, the net effect of the 
reaction is to generate an inhomogeneous temperature gradi- 



ent in the vicinity of the Janus particle, giving rise to propul- 
sion by a thermophoretic mechanism. This scheme is easily 
generalized to allow for different interaction potentials for the 
chemical species on the different hemispheres of the Janus 
particle. In such a situation both thermophoretic and diffusio- 
phoretic mechanisms will operate, and they can act in either 
the same or opposite directions. For the dissociation reaction, 
it is the increased number of product molecules produced in 
reaction, even when the potential energies of interaction of re- 
actant and products with the Janus particle are the same, that 
leads to propulsion by self-diffusiophoresis. Again if these po- 
tential energies were different the propagation velocity would 
change. 

The work presented here provides a method that allows 
one to explore self propulsion for a variety of motor geome- 
tries and catalytic reaction mechanisms at a mesoscopic level. 
Through such studies one can investigate motor dynamics on 
scales where the validity continuum theories should be tested, 
and study situations where analytical methods are difficult to 
carry out due to the complex nature of the reaction dynamics. 
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